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Abstract — In the deregulated electricity market, each 
generating company has to maximize its own profit by 
committing suitable generation schedule termed as profit 
based unit commitment (PBUC). This article proposes a 
Parallel Particle Swarm Optimization (PPSO) solution to the 
PBUC problem. This method has better convergence 
characteristics in obtaining optimum solution. The proposed 
approach uses a cluster of computers performing parallel 
operations in a distributed environment for obtaining the 
PBUC solution. The time complexity and the solution quality 
with respect to the number of processors in the cluster are 
thoroughly tested. The method has been applied to 10 unit 
system and the results show that the proposed PPSO in a 
distributed cluster constantly outperforms the other methods 
which are available in the literature. 

Index Terms — deregulated market, profit based unit 
commitment, particle swarm optimization, distributed 
environment, parallel processing, parallel particle swarm 
optimization. 

I. Introduction 

The GENCOs objective is to maximize the profit and to 
place proper bid in the market. In order to do this generation 
companies need to determine the schedule and operating 
points based on the load and price forecasted. The traditional 
unit commitment problem objective is minimizing the cost of 
operation subject to fulfillment of demand. However in a 
deregulated environment the traditional unit commitment 
objective needs to be changed to maximize the profit with 
relaxation of the demand fulfillment constraint. This unit 
commitment is referred to as profit based unit commitment. 

[1]- 

A competitive and deregulated framework is replacing 
traditional and centralized regulation in many electric power 
systems around the world. With the promotion of deregulation 
of electric power systems, operation, planning and control 
aspects in traditional power system need to be changed [2- 
3]. In this new paradigm, the signal that would enforce a 
unit's on/off status would be the price, including the fuel 
purchase price, energy sale price, ancillary service sale price, 
and so on. There are many solution techniques such as 
integer programming; dynamic programming, Lagrangian 
relaxation and genetic algorithms are available to solve the 
PBUC problem [4-6] . Researchers also presented a review on 
deterministic, meta-heuristic and hybrid approaches of 
generation scheduling in both regulated and deregulated 
power markets [7] . 
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The basic idea of LR is to relax the UCP constraints into a 
small sub-problem, which is much easier to solve, and then 
coordinated by a master problem via properly adjusting a 
factor called Lagrangian multiplier. For all that, it has proven 
to be a very difficult task that may come sometime from oscil- 
lation of their solution by slight change of the multiplier. In 
order to overcome these complex mathematical problems, there 
are other method of computational methodology, which is 
shared by popular artificial intelligence such as genetic algo- 
rithm and evolutionary programming. 

Charles W. Richter et.al presented a PBUC problem 
formulation using genetic algorithm (GA) which considers 
the softer demand constraints and allocates fixed and 
transitional costs to the scheduled hours [6]. Pathom 
Attaviriyanupap et.al proposed a method that helps GENCO 
to make a decision on how much power and reserve that 
should be sold in markets, and how to schedule generators 
in order to receive the maximum profit [7] . Here the authors 
have considered both power and reserve generation at the 
same time. In [8], H. Y. Yamin et.al proposed an auxiliary hybrid 
model using LR and GA to solve UCP. GA is used to update 
the Lagrangian multiplier also presented their view on the 
profit based unit commitment in day- ahead electricity markets 
considering the reserve uncertainty [9]. 

The optimization method known as particle swarm 
optimization (PSO) algorithm developed by Eberhart and 
Kennedy is successfully applied to solve nonlinear 
optimization problems. Therefore an attempt is made to solve 
the PBUC problem using this algorithm. The swarm-based 
algorithm described in this paper is a search algorithm capable 
of locating optimal solutions efficiently. 

The proposed method is applied to solve PBUC problem 
with ten-unit and hundred-unit-test systems. The performance 
of the PPSO algorithm in terms of solution quality is compared 
with that of other algorithms reported in literature for the 
above mentioned problem in power system operation. 
Likewise, simulation results demonstrate the feasibility and 
effectiveness of the proposed method, as compared with the 
results available in the literature. 

II. Problem formulation 

The objective of the PBUC problem is to obtain the optimal 
unit commitment schedule thereby maximizing GENCOs profit. 
The problem formulation is given as follows: 

Maximize PF = RV -TC (1) 

or 



24 



-ACEEE 



ACEEE Int. J. on Electrical and Power Engineering, Vol. 02, No. 03, Nov 201 1 



Minimize TC - RN (2) 

where pp is the total profit ($), rv is the total revenue ($) 
and TC is the total cost ($). 
Here, 

TC =f J it[C l (P {lJ) I (IJ) ) + ST l ] (3) 



(=1 !=1 

T N 



*v=ZZ [(<x,(^ w / (W )] 



(4) 



r=l i=l 



where C- is the production cost which is calculated by using 
the equation (10). P( i>t )te the power level of \ th generator 
unit at t ,h hour (MW), 7 ( , t) is the commitment state of / * unit 
at f' h hour, ST t is the startup cost ($), t is the index for time, 
X is the dispatch period in hours, { is the index for generator 

unit, N is the total number of generating units and a (f) is 
the forecasted market price for energy at time t. 

C i (P ( , l) ) = a + b*P (U) +c*Pl ) (5) 

where a , t> ar, d C are the fuel cost co-efficients. 
System Constraints 
Demand constraints 

i> (W / (i , t) <A r = i....J (6 ) 

/=i 

where D, is the total system demand at time t. 
Here, demand and reserve constraints are different from 
traditional UC problem because GENCO can now select to 
produce demand and reserve less than forecasted level if it 
creates more profit. 
Unit constraints 
1. Unit power limit 

P <P I <P <7"> 

1 i.rain — ' (i,t) * (!*,() — 1 /.max V ' > 



where P inia is the minimum power output of \ th generator 

unit (MW) and P i:mx is the maximum power output of 

generator unit (MW). 

2. Minimum Up and Down time constraints 

[X "" (i,t- 1) -T m (i)] * [/<,,,_„ - /,,,,]> (8) 

[X ° ff (i,t- 1) -T off (/)] * [I (lt) - 7 ( ,,,_ n ] > (9) 

where X°"(i,t) is the "On" duration of /"' generator unit till 

timet, X off (i,t) is the "Off duration of [ ,h generator unit till 

time t, T on (i) is the minimum up-time of f th generator unit 

and T off (i) is the minimum down-time of i' h generator unit. 



III. Parllel clustering environment 

In clusters, powerful low cost workstations and/or PCs 
are linked through fast communication interfaces to achieve 
high performance parallel computing. Workstation clusters 
have become an increasingly popular alternative to traditional 
parallel supercomputers for many workloads requiring high 
performance computing. The use of parallel computing for 
scientific simulations has increased tremendously in the last 
ten years, and parallel implementations of scientific simulation 
codes are now in widespread use [10, 11]. There are two 
dominant parallel hardware/software architectures in use 
today are Distributed memory, and Shared memory. In shared 
memory systems, parallel processing occurs through the use 
of shared data structures, or through emulation of message 
passing semantics in software. Distributed memory systems 
are composed of a number of interconnected computational 
nodes, which do not share memory, but can communicate 
with each other through a high-performance ether net switch 
(HPES) as shown in Figure 1. Parallelism is achieved on 
distributed memory systems with multiple copies of the 
parallel program running on different nodes, sending 
messages to each other to coordinate computations. The 
cluster should perform as a parallel computing resource, 
achieving higher performance than possible using 
workstations configured in a more standard way. The nodes 
in the cluster are always used in groups, not individually as 
in a general purpose workstation laboratory. 

Speedup factor and efficiency: 

To evaluate the parallel performance of the PPSO 
algorithm, the speedup factor SW h and efficiency EW h of the 
cluster [12-13] is calculated as follows; 

SW„=W,/W hl (10) 

EW h =SWjW h (11) 

where W and W h are the execution time of single processor 
and cluster respectively. 
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Figure 1. Distributed cluster of workstations (20 Nodes) 
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III. Particle swarm optimization algorithm 

Particle Swarm Optimization (PSO) is an optimization 
technique inspired from bird flocking which is developed by 
Dr Eberhart and Dr. Kennedy way back in the year 1 995 [ 1 4] . 
It is a population based algorithm where each individual 
(particle) in the population is a potential solution, flies in the 
D dimensional problem space with a velocity which is 
dynamically adjusted according to the flying experiences of 
its own and its colleagues. 

A) Standard PSO algorithm 

Suppose that the search space is D-dimensional and m 
particles form the colony. The z" 1 particle represents a D 
dimensional vector X.(i=l, 2. . . in). It means that the i ,h particle 

positions at X ( =(x n ,X i2 , >■*#)) (i=l, 2... m) in the 

searching space. The position of each particle is a potential 
result. The calculation of the particle's fitness is carried out 
by putting its position value into a designated objective 
function. When the fitness is higher, the corresponding X is 
"better". The z" 1 particle's "flying" velocity is also a D- 

dimensional vector, denoted as V i = (v n ,V i2 , j v ,£)) 

Denote the best position of the z' th particle as 

P. = (p n , p j2 , 'Pid^ anc ^ trie b est position of the 

colony as P g =(p gl , P gl , >P s d) respectively. The 

PSO algorithm can be performed by the following equations 
(12,13). 

X d (k+l)=X d (k>c 1 r 1 (P Id (k>x ld (k)>c 2 r 2 (P gd (k>x Id (k) (12) 

x id (k + \) = x id {k) + v id (k+\) (13) 

Where k represents the iterative number, 
c 1? c, are learning factors. Usually c [ =c,=2 r , r, are random 
numbers between (0, l).The termination criterion for the 
iterations are determined according to whether the max 
generation or a designated value of the fitness of P is reached 

B) Parallel PSO algorithm 

PPSO algorithm is implemented to determine the 
commitment status of each unit over a scheduled period of 
(24 hours) time in order to maximize the profit. The procedure 
of the proposed algorithm to solve PBUCP is as follows. 
Step 1 : Generator and PSO Parameters Specification 

Specify the generator minimum and maximum generation 
limits, minimum up and down time constraints and start up 
cost of each unit. Specify the PPSO parameters such as 
population size (M), inertia weight factor (w), dimension of 
the system (D), acceleration constants (c { and c,), velocity 
maximum and minimum limits, maximum iterations (Max iter). 
Set iteration number iter=l and time t=l . 
Step 2: Particles sharing policy 

Master node decides the sharing of particles by particles 
sharing policy (PSP). Therefore the number of particles 
allocated in a slave processors or workers is given by 



Ix, if particles divisible by W/ t 
x + lforthe first hx slaves and 
x for the remaining , Otherwise 

where 

x= floor {N p artic i es IW h ) 
hx=N Par1kles -(x*W h ) 

Master node allocate (x+1) particles to the first hx slaves in 
the Wh cluster (Wi..Wf lx ...Wf 1 ) and x ants to the remaining 
slaves (W hx+l ...W h ) . Where JV pmi . rfa =population. 

("start 



(14) 

(15) 
(16) 



J 



Inifelize generator and PPSO paraa ater paraa e ters 

such as inerte weight fac tor. acceleration factors* 

mai and mill velocity fails, population size and 

no bar of workers, iter=Q. 



Mesfer shares tie partibs a-nong the 
w offers processors by PSP 




Genera fe initiate solution for each individual in the 
s^ara us in; uniforaty distribufed random function 



Repair each individual for a mtri ua> up and down 
true constraint vctstbr. 



Coapute profit foi each individual in the swa 
using objx rive f unc liar. 



Cocd pare each particle profit with its own tea t 
position (Pbesf). if new position is be tfer than the 
old one. set the new position as PbesL otherwise 

retain old one as PbesL The best Pbesl(BP) is 
sends to the m aster. 




Master compares the best Pbest re:eived froa the slaves and 
sebc t the bes t global best (GP) so far obtained and requests 
= " -■■ a-r^iers :: i: zizqsl^pistr.g 




Figure 2. Flowchart of PPSO for PBUC 

Step 3: Initialization of Individual in the swarm 

The initial solution of each individual U=[u , u , ... u J, (j=l, 

j L nl n2 nT J ^ ' 

2, ...,M), (n=l, 2, , N)) for complete M population is 

generated randomly. The position of each unit u of each 
particle is generated using uniformly distributed random 
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function, which generates either or 1 . Similarly the initial 
velocity of each particle is generated randomly using 
uniformly distributed random function, which generates a 
real value between V and V .The representation of each 

nun max x 

individual for 'N' number of generating units for a scheduled 
period of time is as follows: 



u = 



u. 



u 



21 



'12 



"22 



U 



T\ 



T2 



"IN 



'IN 



"TN 



(17) 



Step 4: Defining the evaluation function 

The merit of each individual particle in the swarm is found 
using a fitness function called evaluation function. Each 
particle in the population is evaluated using the objective 
function given by ( 1 ). 

Step 5: Repair Minimum up and down time constraints 
violation 

Repair each unit for each particle in the swarm for minimum 
up and down time constraints violation. 
Step 6: Modifying best particle position (Pbest) 

To modify the position of each individual in the next stage 
is obtained from equation (12). 
The weighting function is defined as follows 



iter m „ 



iter 



(18) 



Where, 



iter 



Initial, final weights 



- Maximum iteration number 

iter - Current iteration number 

To control excessive roaming of particles, velocity is restricted 

between - v" 1 "" and +V max . 

The maximum velocity limit for the j th generating unit is 

computed as follows: 

p max _ p min 

v m «=ii J. (19) 



The particle position vector is updated using equation (13). 
The values of the evaluation function are calculated for the 
updated positions of the particles. Evaluate each particle using 
object function. Compare each particle evolution value with 
its own best position (Pbest or BS). If the present particle 
position is better than the old value set new particle position 
as Pbest, otherwise retain old value. 
Step 7: Computation of global best 

Master receives the information of local best solutions 

(BS r BS 2 , ,BS h ) from the workers and computes the best 

solution among them as the global best. Whenever a global 
best solution is selected by the master, and if the total profit 
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is found to be more than the maximum total profit computed 
so far, then the present global best is memorized, or else the 
previous maximum total profit solution is retained as global 
best. The new global best is sent to the workers and the 
workers saved the received global best as their global best. 
Step 8: Memorize the best solution obtained so far and 
increment iteration number. Stop the process if iteration 
number is equal to the maximum number of cycles. Otherwise 
go to step 4. 

Step 9: Increment the time and repeat step 3 to step 8 for the 
given scheduled period (24 hours) of time. 
The flowchart of the proposed method is shown in the Figure 
2. 

IV. Numerical results 

The PABC method for PBUC is first tested on 10 unit 
system available in the literature as Case 1 . It is also validated 
on multiple test systems of 100 units in Case 2. The parallel 
computation is carried out in the MATLAB® environment of 
R2007b using distributing computing toolbox. The parallel 
computation is carried out through distributed memory 
environment. In a distributed environment, a cluster with the 
maximum size of 20 nodes/processors (Pentium - IV 3.40 GHz, 
1GB RAM) is used. 

A. Case 1: 10 Unit System 

In order to participate in the market, GENCOs have to 
prepare a self commitment according to the forecasted load 
and price. In this case, the commitment schedule is prepared 
to maximize the GENCOs' profit by calculating the generator 
coefficients with the satisfaction of constraints. Here the profit 
of the company gets the first priority and the demand 
satisfaction is not mandatory. So, GENCOs will make the self 
commitment depending upon the forecasted price to get 
surplus profit. The test system consists of 10 generating 
units. Here the generating unit data and load data are taken 
from [15]. The constraints included for PBUC in [15] are 
considered. Based on the forecasted market price of energy 
information, the proposed PPSO model is used to generate 
dispatch schedule for 24 hours time horizon. The dispatch 
schedule of ten unit system is given in TABLE I. Optimal 
parameters obtained by trial and error for PPSO is as follows: 
Population size=200, Acceleration coefficients, c 1=0.2, c2=0.2, 
Inertia weight: W =0.9, W =0.2 and Maximum 

max min 

iteration=300. 

The comparison of the proposed method with other 
existing methods given in TABLE II proves that PPSO gives 
better solution, i.e. a difference in profit of $1260.23 is 
achieved when compared to Muller method [15]. The time 
taken to get the best schedule is 168 sec. The PPSO yields a 
higher profit of 1 .2 % than the Muller method. Figure 3 shows 
the execution time achieved by the different cluster sizes. 
The execution time will reduce, when the cluster size 
increases. TABLE III shows the speedup factor and efficiency 
achieved by different sizes of cluster. When the cluster size 
increases the speedup factor also increases, i.e., performance 
of the cluster will increase, when its size increases. 
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Table i. Dispatch schedule for 10 unit system 



Hr 


Load 
(MW) 


in 

(MW) 


U2 
(MW) 


U3 

(MW) 


U4 
(MW) 


U5 
(MIS 7 ) 


U6 
[MW) 


U7 

(MW) 


us 

(MW) 


U9 
(WW 7 ) 


U10 
(MW) 


Fuel Cost 

($) 


Transition 
Cost (S) 


Total Com 


Revenue 
(5) 


Profit ($) 


1 


700 


455 


245 


























13683.13 





13683.13 


15505.00 


1S21.S7 


7 


750 


455 


295 


























14554.50 





14554.50 


16500.00 


1945.50 


3 


8:0 


4:5 


395 


























16301. S9 





16301. S9 


19635.00 


3333.11 


4 


950 


455 


455 


























17353.30 





17353.30 


20611.50 


3258.20 


5 


1000 


455 


455 


























17353.30 





17353.30 


21157.50 


3S04JO 


6 


1100 


455 


3S5 


130 


130 




















21S79.33 


2220 


24099.33 


25245.00 


1145.67 


7 


1 ISO 


455 


435 


130 


130 




















22755.04 





22755.04 


25S75.00 


3119.96 


S 


1200 


455 


455 


130 


130 




















23105.76 





23105.76 


25915.50 


2S09.74 


9 


1300 


455 


455 


130 


130 


130 

















26184.02 


1800 


27984.02 


29640.00 


1655.98 


10 


1400 


455 


455 


130 


130 


162 


6S 














2S76S.21 


340 


2910S.21 


41090.00 


119S1.79 


11 


1450 


455 


455 


130 


130 


162 


80 














29047.98 





29047.9S 


42571. SO 


13523.S2 


i: 


1500 


455 


455 


130 


130 


162 


80 














29047.98 





29047.9S 


446S9.S0 


15641.32 


13 


1400 


455 


455 


130 


130 


162 


6S 














2S76SJ1 





2S76S.21 


34440.00 


5671.79 


14 


1300 


455 


455 


130 


130 


110 


20 














26588.96 





265SS.96 


31S50.00 


5261.04 


15 


1200 


455 


455 


130 


130 




















23105.76 





23105.76 


26325.00 


3219.24 


16 


1050 


455 


335 


130 


130 




















21005.17 





21005.17 


13415.00 


2409.83 


17 


1000 


455 


235 


130 


130 




















20132.56 





20132.56 


22250.00 


2117.44 


IS 


1100 


455 


3S5 


130 


130 




















21S79.33 





21879.33 


24255.00 


2375.67 


19 


1200 


455 


455 


130 


130 




















23105.76 





23105.76 


25974.00 


2S6SJ4 


20 


1400 


455 


455 


130 


130 




















23105.76 





23105.76 


26500.50 


33 94.74 


21 


1300 


455 


455 


130 


130 




















23105.76 





23105.76 


27027.00 


3921.: 4 


22 


1100 


455 


385 


130 


130 




















21879.33 





21879.33 


25245.00 


3365.67 


23 


900 


455 


445 


























17177.91 





17177.91 


20475.00 


3297.09 


24 


SOO 


455 


345 


























15427.42 





15427.42 


1S04Q.00 


2612.58 


Total (S) 


525316.37 


4360 


52967637 


634232.60 


104556.23 



Table ii. Comparison of pbuc solutions (10 unit system) 



Method 


Profit (S) 


TS-RP[16] 


1010S6.00 


TS-IR.P[16] 


10326100 


Muller methodfl 5] 


103296.00 


PSO 


104356.00 


PPSO 


104556.23 



J SCi 




B. Case 2: 100 Unit Syatem 

This test system consists of multiple generating units 
such as 100 generating units. More number of generating 
units is considered in order to validate the feasibility of the 
application of PPSO for large scale power system. The data 
for different groups of generating units are obtained by 
duplicating the 10 unit system data. The demand is multiplied 
with respect to the system size; however the generating limits, 
the minimum up/down time constraints remain the same. 
Based on the forecast market price of energy information, the 
proposed PPSO model is used to generate dispatch schedule 
for 24 hours time horizon. The parameter setting of the 10 
unit system is extended for the multiple test systems. 



::: 



2 4 6 S 10 12 14 16 1B 20 

Number of Processors 

Figure 3. Execution time chart for 10 unit system 

Table in. Comparision of speedup factor and cluster efficiency (10 unit 

system) 




w k 


2 


4 


6 


S 


10 


12 


14 


16 


IS 


20 


SpeeAp 
factor 


1.6 


2.69 


4.01 


5.04 


6.03 


7.01 


7.49 


8.11 


S.43 


S.9S 


EWtfK) 


S0.01 


6725 


66.S3 


63.01 


60.03 


5S.41 


53.51 


50.6S 


46.S3 


44.91 



6 S 10 12 14 
NuberofPro cesso rs 
Figure 4. Execution time chart for 100 unit system 
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Table iv. Commitment status for 100 unit system 



Hr 


Unit &tatu& 


Fuel Cost 

($) 


Tramiti 
on Co&t 


Total 
Cost(S) 


Revenu 
e{$) 


Profit (S) 


1 


nooooooooi 100000000110000000011000000001100000000110000 do do 

1 1 00000 0001 1 000 00 000 1 1 000 00 000 1 1 00000000 


136S313 







136831.3 



155050 


18218.70 


2 


1 1 000000001 1 00000000 1 1 00000000 1 1 00000000 1 100000000 1 100000000 
1100000000110000000011000000001100000000 


145545.0 







145545.0 



165000 


19455.00 


3 


110000000011000000001100000000110000000011000000001100000000 
1100000000110000000011000000001100000000 


163018.9 






163018.9 



196350 


33331.10 


4 


11 00000 0001 100000000 11 00000000 11 000000 00 11000 000 00 1100O0 0000 
1 1 00 00 1 1 000 00 000 1 1 000 00 000 1 1 000 000 00 


173533.0 






173533.0 



206115 


32582.00 


5 


11 00000 0001 100000000 11 00000000 11 000000 00 11000 000 00 110000 00 00 
11 00000 0001 100000000 11 00000000 11 00000000 


173533.0 







173533.0 



211575 


38042.00 


6 


1 1110000001111000 000 11 11000 00011110000 0011110000001111000000 
1 1110000001111000 00011 11000 00011 110000 00 


218793.3 

1 


22200 


240993.3 



252450 


11456.69 


7 


111100000011110000001111000000111100000011110000001111000000 
111100000011110000001111000 00011110000 00 


227550.4 

1 





227550.4 



258750 


31199.59 


£ 


111100000011110000001111000000111100000011110000001111000000 
1111000000111100000011110000001111000000 


231057.5 
9 





231057.6 



259155 


28097.41 


9 


111110000011111000001111100000111110000011111000001111100000 
1111100000111110000011111000001111100000 


261S40.2 

1 


1S0OO 


279840.2 



296400 


16559.79 


10 


111111000011111100001111110000111111000011111100001111110000 
1111110000111111000011111100001111110000 


287682.1 
3 


3400 


2910S2.1 



410900 


119817.9 



11 


111111000011111100001111110000111111000011111100001111110000 
1111110000111111000011111100001111110000 


290479.7 

s 





29Q479.S 



425718 


13 523 82 




12 


111111000011111100001111110000111111000011111100001111110000 
1111110000111111000011111100001111110000 


290479.7 
8 





2904 7 9. S 



446898 


156418.2 



13 


111111000011111100001111110000111111000011111100001111110000 
1111110000111111000011111100001111110000 


287682.1 

3 





2S76S2.1 



344400 


56717.87 


14 


111111000011111100001111110000111110000011111000001111100000 
1111100 0001111100 00011 11100 00011 11100000 


263368.0 






26336S.0 



31S5O0 


55132.00 


15 


1 1110000001111000 000 11 11000 00011110000 001111000000111100 00 00 
11110000001111000 0001111000 00011 11000000 


231057.5 
9 





231057.6 



263250 


32192.41 


16 


111100000011110000001111000000111100000011110000001111000000 
1 1110000001111000 00011 11000 00011 11000000 


210051.7 

1 





210051.7 



234150 


2409829 


17 


1111000 0001111000 0001111000 00011110000 001111000000111100 00 00 
1111000 0001111000 000 11 11000 00011110000 00 


201325.6 

1 





201325.6 



222500 


2117439 


IS 


1111000 0001111000 000 11 11000 00011110000 0011110000 001111000000 
1111000000111100000011110000001111000000 


218793.3 

1 





2187933 



242550 


23756.69 


19 


111100000011110000001111000000111100000011110000001111000000 
1111000000111100000011110000001111000000 


231057.5 
9 





231057.6 



259740 


28682.41 


20 


1 1 1 1 0000001 1 1 1000000 1 1 1 1000000 1 1 1 1000000 1 1 1 1000000 1 1 1 100O0O0 
1111000 0001111000 0001111000 0001111000000 


231057.5 
9 





231057.6 



265005 


33947.41 


21 


111100000011110000001111000000111100000011110000001111000000 
1111000000111100000011110000001111000000 


231057.5 
9 





231057.6 



270270 


39212.41 


22 


111100000011110000001111000000111100000011110000001111000000 
1 1110000001111000 00011 11000 00011 110000 00 


218793.3 

1 





218793.3 



252450 


33656.69 


23 


11 00000 0001 1000 00 000 11 000 00 000 11 000 000 00 1100000000 1100000000 
1 1 000000001 1 00000000 1 1 00000000 1 1 00000000 


171779.1 







171779.1 



204750 


32970.90 


24 


11 00000 0001 1000 00 000 11 000 00 000 11 000 000 00 1100000000 1100000000 
1 1 000000001 1 00000000 1 1 00000000 1 1 00000000 


154274.2 







1542 74 2 



1S0400 


26125.80 


Total (5) 


5250642. 


43600 


5294242. 


634232 


104S0S3. 



Table v. Comparision of speedup factor and cluster efficiency (100 unit 
system) 



w„ 


2 


4 


6 


g 


10 


12 


14 


16 


IS 


20 


Speedup 
factor 


1.S5 


3.01 


423 


5.71 


6.9S 


S.14 


S.62 


928 


10.39 


11.05 


EW a (%) 


92.51 


7525 


70.51 


7137 


69 Jl 


67. S3 


61.57 


5S.01 


57.72 


5525 



The Commitment status of 100 unit system is given in TABLE 
IV. Figure 4 shows the execution time achieved by the differ- 
ent cluster sizes. The execution time will reduce, when the 
cluster size increases. TABLE V shows the speedup factor 

©2011 ACEEE 
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and efficiency achieved by different sizes of cluster. When 
the cluster size increases the speedup factor also increases, 
i.e., performance of the cluster will increase, when its size 
increases. When a single processor is used, it consumes 
more time for execution, i.e., 1728.01 sec for 100 units. The 
execution time for the 20 node cluster of 100 units are around 
156.38 sec. It clearly shows the execution time decreases as 
the number of processor increases. Each test system has 
been tested for 30 trial runs and the best results are pre- 
sented. 
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V. Conclusions 

This paper proposes a MPI based PPSO model for PBUC, 
computing in parallel, in a distributed environment. The 
approach is simple, efficient, and economic and can be 
extended for making smarter decisions in a large scale power 
system. Simulation results obtained from the cluster 
demonstrate the accuracy of the proposed algorithm and its 
capability of greatly reducing the execution time. 
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